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Abstract 

Measurements of the lowest-lying eigenvalues of the staggered fermion Dirac operator are made 
on ensembles of equilibrium gauge field configurations in quenched SU(3) lattice gauge theory. The 
results are compared with exact analytical predictions in the microscopic finite- volume scaling regime. 
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The study of the eigenvalue spectrum of the Dirac operator in QCD and related gauge theories has a long 
history, most notably due to the relation between zero modes, near-zero modes, topology, and chiral 
symmetry breaking. More recently the subject has undergone a strong revival. This was originally 
sparked by the observation made by Leutwyler and Smilga [Q] that in sectors of fixed topological 
charge v one can derive exact analytical results for the distribution of small eigenvalues of the Dirac 
operator. The idea is that there exists a finite-volume scaling regime for the euclidean gauge theory 
where the partition function can be computed exactly. One crucial input is the assumption that chiral 
symmetry is spontaneously broken according to the conventional scenario. The finite-size scaling regime 
is unphysical in the sense that it restricts the associated pseudo-Goldstone bosons to be very light, so 
light that only their zero-momentum modes contribute significantly to the euclidean path integral. The 
conditions for this simplification are two- fold. First, the euclidean finite- volume partition function of 
QCD must be accurately described by the effective hadronic excitations, rather than the underlying 
quark and gluon degrees of freedom. This imposes the large- volume condition 1/Aqcpj <C L, where 
L is the linear extent of the four-volume. Second, this linear extent L should always be insignificant 
compared with the Compton wavelength of the pseudo-Goldstone bosons: L <C While this latter 

condition is inappropriate for actual, physical, no-so-light pions, it is of course just the right limit for 
studying spontaneous chiral symmetry breaking. More importantly in this context: it is ideally suited 
for finite-volume studies of lattice gauge theory. 



A major breakthrough in the understanding of the Dirac operator spectrum in the above finite- volume 
region has been the suggestion by Verbaarschot and co-workers that random matrix theory can give 
exact analytical predictions for the spectrum in a suitably defined microscopic regime ||, [|. Gauge 
theories based on various gauge groups and with fermions in either fundamental or adjoint representa- 
tions fall into essentially just three different universality classes, which beautifully fit into corresponding 
classifications of random matrix theory ensembles Q. In this paper we shall restrict ourselves to SU(3) 
gauge theory with quarks in the fundamental representation of the gauge group, where universality of 
all microscopic correlation functions, including the microscopic spectral density itself, has been proven 
|| . Moreover, it has recently been shown that all these exact results can also be derived directly from 
the finite- volume partition function of QCD |Jsj, ||] , without recourse to random matrix theory. Also the 



universal smallest-eigenvalue distribution can be computed in this way [1C] 



So far the analytical predictions for the microscopic Dirac operator spectrum have been tested in 4- 



dimensional SU(2) lattice gauge theory with staggered fermions [11|, and in 3-dimensional SU(3) lattice 



gauge theory with staggered fermions [12|. There have also been preliminary studies of the microscopic 
spectrum in 2-dimensional U(l) lattice gauge theory with a fixed-point Dirac operator |]13|| . However, 
none of these cases test the most important universality class of all, namely that of (in the language 
of random matrix theory) the chiral Unitary Ensemble. This is the universality class of 4-dimensional 
SU(N C ) gauge theories with N c > 3, and quarks in the fundamental representation. 

In this paper we shall present the first results for a lattice gauge theory determination of the microscopic 
Dirac operator spectrum of SU(3) gauge theory with a staggered fermion Dirac operator. The biggest 
advantage of using a staggered Dirac operator is that the microscopic spectrum at any value of the 
gauge coupling should coincide with the continuum predictions, since the universality class of this 
lattice-regularized Dirac operator in this case coincides with the universality class of the continuum 
Dirac operator (for fermions in the same color representation). The main disadvantage is related to the 
fact that all comparisons with analytical results should be done in sectors of fixed topological charge v. 
It is well-known that although a remnant of chiral symmetry is exactly preserved in the lattice theory 
with staggered fermions, the interplay between gauge-field topology and fermions is a subtle issue in 
this formulation. In particular, the existence of v fermion zero modes in a gauge field configuration 
with topological charge v is jeopardized. Earlier studies of the microscopic Dirac operator spectrum 
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in SU(2) lattice gauge theory with staggered fermions [11] have suggested that to a large degree the 
effects of non-trivial gauge field topology are insignificant at feasible values of the gauge field coupling, 
and hence can be ignored at this stage. While we feel that this issue still deserves further clarification, 



we shall here take the same stand as in ref. [11], and simply sum over all gauge field configurations, 



irrespective of topology. Comparison should, in the same spirit, be made only with analytical v = 
predictions. 

Before considering the numerical algorithm that computes the Dirac operator eigenvalues in a given 
gauge field configuration, it is instructive to estimate the number of small eigenvalues it is meaningful 
to include in our study. The crucial observation here is that in the microscopic scaling region near 
A ~ 0, the Dirac operator spectrum will have an average level spacing of order unity. Specifically, we 
shall be concerned with the rescaled, microscopic, spectral density 

where £ = = ?r/9(0) , and p(A) denotes the usual (macroscopic) spectral density. We choose a 

convention where the chiral condensate is defined, with a summation over the number of colors N c , by 
(for one flavor of mass m) 

- viL z(m> • < 2) 

Here Z{m) is the finite-volume partition function. The macroscopic spectral density is normalized so 
that p(X)d\ is the mean number of eigenvalues per unit volume, as in ref. ||]. The exact, universal, 
analytical prediction from both random matrix theory || and the finite-size partition functions ||, || 
reads 

Ps(C) = np{0)^ J Nf +u{() 2 ~ JNf+v+l(z)J Nf+v -i{C,) , (3) 

with C = AT^S = irp(0)V\. This analytical expression has the following qualitative features: It vanishes 
at C, = (consistent with the fact that at any finite volume there is no spontaneous chiral symmetry 
breaking), rises sharply, and then undergoes exponentially damped oscillations around the constant 
value p(0)[]. The oscillations take place on a scale of order A£ ~ 1, and, in fact, correspond directly 
to fluctuations over distinct sharp peaks of the first few individual eigenvalues. While the analytical 
prediction is exact (can be achieved to any degree of accuracy) in the infinite-volume limit with L <C 
l/m-v, it is clear that at any finite volume the prediction will fail eventually. In any case, since the 
distinct feature of the prediction, the oscillations, die out exponentially, we are in fact only interested 
in averaging over the first very few eigenvalues. It would thus be considerable overkill to compute, for 
each gauge field configuration, more than a few lowest-lying eigenvalues. We use the Ritz functional 
method [14] to compute the lowest 10 eigenvalues (and eigenvectors) of — If) 2 , with If) the staggered 



Dirac operator 

P*,v = o H ( u n( x ) s *+w - Uj t (y)6 x ,y +t ^ = Ip efi + p , e . (4) 

p- 

Here f]^{x) = (— l)^"^ x " are the usual staggered phase factors. Introducing e(x) = (— 1)^"> X " , we 
have indicated that If) connects even, e(x) = 1, with odd, e{x) = —1, sites and vice versa. 



For clarity, we recall some well known properties of the staggered Dirac operator here, see e.g. ref. | |15| . 
The staggered Dirac operator is antihermitian and therefore —If) 2 is hermitian and positive (semi-) 
definite. e(x) plays the role that 75 plays for the continuum Dirac operator, e.g. {p,e} = 0, and it 
is easy to see that the eigenvalues of If) come in pairs ±iA. If ip(x) is the eigenvector with eigenvalue 

lr The microscopic spectral density can be viewed as a blow-up of the point p(0) due to the finite volume. We should 
thus have lim<;_ 00 p s (C) = p(0), a matching condition which indeed is satisfied. 
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5.1 
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5.55 


6 


4000 


0.0003 


0.192(2) 



Table 1: Our estimates of p(0) for our various gauge field ensembles. The third column gives the 
number of configurations analyzed, and the forth column the bin size in the histograms. 

iX then e(x)ip(x) is the eigenvector with eigenvalue — iX. It is well known that —If) 2 only couples 
even with even and odd with odd sites. Furthermore, if i\) e is a normalized eigenvector of —Ij) 2 with 
eigenvalue A 2 , non-zero only on even sites, then ijj = jlft 0je ip e is also a normalized eigenvector of —If) 2 
with eigenvalue A 2 , non-zero only on odd sites.0 The eigenvectors of fl with eigenvalues ±iA are then 
i>±{x) = ^ T i^o(x)). 

In any case, it is clear that to obtain the lowest 10 positive (imaginary) eigenvalues of Jp it is sufficient 
to compute the lowest 10 eigenvalues of —]p 2 e , working only on the even sublattice, and then taking 
the square root. This is the approach we have taken here, for which the Ritz functional method is very 
well suited. 

We have generated configurations that are essentially statistically independent by employing a pure 
gauge algorithm of typically 3 to 4 micro canonical overrelaxation sweeps, followed by one heatbath 
update sweep, and repeating this process 10 times between eigenvalue computations. To convert our 
eigenvalue measurements, or actually a histogram of all measured eigenvalues, into a microscopic spec- 
tral density, see eq. (||), we need an estimate of E, or of p(0), in the infinite volume limit. On a finite 
lattice p(X), obtained by normalizing a histogram of the lowest eigenvalues appropriately, will oscillated 
about the infinite volume p(0). We thus obtain p(0) just by averaging the finite volume p(X) over 
the first few (leaving out the very first) oscillations. Our results, together with information about the 
number of measurements and the lattice size, is given in Table @. 

Our results here are restricted to the quenched approximation, which, in the analytical context can be 
achieved simply by taking the Nf — » limit of the formula (H|). While the quenched approximation 
obviously saves an enormous amount of computer time, it also has another unique advantage: The 
massless limit of the Dirac operator is taken automatically. This means that the "valence Goldstone 
boson" in principle, were there no lattice artefacts, would be strictly massless. One boundary of the 
scaling regime, L <C l/m n is thus trivially satisfied from the very beginning. On the lattice, the other 
boundary, 1/Aqcd <C L, is effectively replaced by the condition that the volume should be large in 
physical units. This means that the lattice volume in lattice units must increase steeply with the gauge 
coupling f3 = Q/g 2 if we wish to go towards weak coupling, and continuum scaling. Another way of saying 
that is as follows: to obtain reasonable estimates of p(0) - i.e. for p(A) to be approximately constant 
(apart from the predicted oscillations), we need sufficiently large lattices in physical units. Empirically 
we find that L>2A r TiC (/3) is needed, with A r T)C (/3) denoting the approximate temporal extension of a 
system to be at the deconfinement transition point at given gauge coupling (5. For example, we found a 
6 4 lattice at /3 = 5.7 (iV T)C (5.7) r* 4) and an 8 4 lattice at /3 = 5.85 (iV T)C (5.85) « 6) not to be sufficiently 

2 Due to lattice artefacts the staggered Dirac operator does not have exact zero-modes on equilibrium gauge configura- 
tions at non-zero gauge coupling, g 2 . 
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large. 



We begin our comparison between theory and Monte Carlo results with the distribution of the smallest 
eigenvalue. Here the analytical prediction, which again has been proven to be universal in the random 



matrix theory context, and derivable directly from the finite- volume partition functions too [1C], reads 



^min(C) = vrp(0)|exp[-C 2 /4] • (5) 

Results for different lattice sizes of volume L 4 are shown in the histograms of Fig. ||, which are compared 
with the analytical prediction (||). For each f3- value we have first determined p(0) as described earlier, 
which means that these plots are all parameter-free. The agreement is seen to be excellent indeed. 

Next, we have extracted the spectral density p(\) near A ~ by binning eigenvalues in small intervals 
(see Table |]), an d normalizing the resulting density as described previously. By rescaling the small 
eigenvalues according to £ = irp(0)V\, we obtain the microscopic spectral density as shown in Fig. |2|. 
Again we find absolutely excellent agreement with the exact analytical prediction (|3|) , which is shown as 
the full line. We have the poorest statistics for (3 = 5.1, L = 8 (only 1500 independent configurations), 
and indeed this is the plot for which we get the least impressive agreement. However, taken together, 
over several different /3-values and several different volume sizes, the precise prediction is very well 
reproduced. 



We conclude by noting that the exact analytical predictions for the microscopic Dirac operator spec- 
trum, and the smallest Dirac eigenvalue distribution, now have been successfully confirmed by direct 
numerical computation in SU(3) lattice gauge theory with staggered fermions. The relevant universality 
class is, in the language of random matrix theory, that of the chiral Unitary Ensemble. For gauge group 
SU(N C ) with N c > 3 and fermions in the fundamental representation, it coincides with the universality 
class of the continuum Dirac operator. 
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Figure 1: The distribution of the lowest eigenvalue for our different lattice ensembles of Table [E]. 
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Figure 2: The microscopic spectral density for our different lattice ensembles of Table [l]. 
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